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ANALYSIS OF THE SURFACE LOAD AND RADIATED SOUND 
OF A 'VIBRATING AIRFOIL WITH APPLICATION TO THE 
EXPERIMENT OF BROOKS 


John E. Yates 
SUMMARY 


A method is developed for calculating the surface load and radiated 
sound from a vibrating surface in a compressible viscous fluid. The method 
is applied to a thin two-dimensional elliptic cross-section. For large 
values of the viscous diffusion parameter, o, the surface load tends to an 
elliptic distribution in agreement with the results of inviscid theory when 
edge pressure continuity (the Neumann condition) is enforced. For thin 
surfaces, the surface load is insensitive to variations in the thickness 
ratio. 


A three-dimensional spectral technique is developed to calculate the 
inviscid surface load and radiated sound from a thin vibrating airfoil. 
Calculations are compared with an experiment of Brooks (1977). The 
inviscid theory predicts the correct form of the far field sound pressure 
and its phase. The actual levels are somewhat sensitive to the choice of 
theoretical spanwise surface pressure mode but are in better agreement with 
experiment than the surface pressure. The comparison of theoretical and 
experimental surface pressure indicates that the viscous theory, used to 
validate the inviscid theory, is either inadequate or there is a source of 
experimental error. A previous explanation (Yates, 1980) for the 
experimental data, that there is an equipartition of surface work into 
sound and viscous dissipation, is verified if flat surface loading is taken 
to occur near the airfoil edges. However present calculations fail to 
confirm the type loading indicated' by the experiment. 



1 . INTRODUCTION 


Brooks (Ref. 1) published the results of an experiment that was 
carefully designed to evaluate and hopefully validate the Kirchhoff 
integral formulation of sound radiation from a vibrating airfoil. The 
surface motion and pressure were measured and used in the Kirchhoff 
integral to calculate the far field sound. Simultaneously the sound field 
was measured and compared with the calculated field. The outcome was 
unexpected. The calculated sound field is 3 to 5 dB more intense than the 
measured field (see Appendix A, Figs. A. 4, A. 5). Also, the measured 
surface pressure distribution appears to have a virtual discontinuity 
across the leading and trailing edges although the spanwise distribution 
does decay smoothly to zero at the edges as may be expected (see Appendix 
A, Figs. A. 2, A. 3)* 

Yates (Refs. 2 and 3) published a plausible explanation of Brooks' 
results. Based' on the premise that the measured surface pressure 
distribution is correct, a calculation was made of the partitioning of the 
surface work into radiated sound energy and vorticity formation near the 
edges. At the Helmholtz numbers of the experiment, the surface work is 
almost equally partitioned into sound and vortical energy produced near the 
leading and trailing edges. On the other hand, the inviscid Kirchhoff 
formulation, used by Brooks to calculate the far field, assigns all of the 
surface work to the far field. Hence, the 3 to 5 dB overestimation of the 
sound. Based on the energy argument of Yates, two conclusions were drawn. 

1) Viscous forces cannot be neglected in the Integral relation (e.g., see 
Ffowcs Williams and Hawkings, Ref. 4) between the surface and far field. 

2) Viscous stresses must be included in the calculation of the surface 
load. 

In a subsequent study the present author focused on the aerodynamic 
problem of calculating the surface load. It was shown (Ref. 5) that the 
inviscid aerodynamic problem is not well posed in that it admits 
eigensolutions . Furthermore, the far field sound depends on the magnitude 
of the eigensolution and is also not unique. By relaxing the usual Neumann 
edge uniqueness condition (continuity of pressure), it was possible to find 
by arbitrary adjustments of the eigensolution a surface pressure 
distribution in much better agreement with the one measured by Brooks. 
Because these inviscid results were encouraging, a direct viscous 
calculation of the load on a vibra ting fla t plate was made and included in 
Ref. 5. For Stokes number* (o - /wc V4v) of order 10, a singular edge 
loading is calculated. For o * 100, the edge load decays smoothly to zero 
in accordance with the Neumann condition. The Stokes number of the Brooks 
experiment is approximately 3000 (see Appendix A), and the calculated 


*See remarks on P. 11 concerning this terminology. 
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viscous edge loading is in complete agreement with the inviscid theory plus 
Neumann condition. Based on the results of viscous thin airfoil theory 
(Ref. 6), Yates suggested that the viscous calculations must be made for 
the actual cross-sectional geometry of the airfoil with Stokes number based 
on the local edge radius. 

The purpose of the present study is to perform the suggested detailed 
calculations for the surface pressure distribution and also for the sound 
field. Based on the results, evaluations are made of the experiment and 
its previous explanations. It is shown here that even with finite leading 
and trailing edge radii, the viscous theory predicts a surface load 
distribution that is in agreement with results of inviscid theory plus the 
Neumann edge condition. This is taken to indicate that either the viscous 
theory as applied is inadequate to predict the experimental results, or 
something was overlooked in the experiment. The present report does not 
offer a conclusive answer to this alternative. However, to conclude our 
study, specific questions were examined in detail with the use of inviscid 
load and acoustic theory. Follow up experimental and theoretical efforts 
are proposed in order to resolve unanswered questions of fundamental 
importance in our understanding of the role of viscous stresses in surface 
loading and noise production. 
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NOMENCLATURE 


a 

A 


A* 

b 

c 


c 


P’ 


c 


V 


C mn 

d n (s), d t (s) 
D n (s), D t (s) 

%» d T 
e, e' 

E(k), K(k) 

F n 

F(w) 


background speed of sound 

b/c, aspect ratio for rectangular airfoil 

coefficients in the series solution of the 
thickness problem see (2.3-16) 

effective, aspect ratio, see Section 3- 1 * 

airfoil span, see Fig. 3-1-1 

airfoil chord 

specific heats at constant pressure and 
constant volume 

matrix coefficient, see (2.3- 18) , (2.3.21 ) 
or (3-2. 17), (3-2. 23) 

displacement vector, see (2.2.1) 

see (2.2.8) 

see (2.2.18) 

see (2.3.4) 

see (3.5.1), (3-5.3), (3-6.1). (3-6.2) 

standard elliptic integrals, see (2.3-28) 

see (2.3-24), through (2.3-27) 

conformal mapping function, see (2.2.19) 
and (2.3.1) 
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Bf 8 0 


8m 

G m 

h o 

h* 

fi(x) 


H 

H f 

H„ 


i, j, i< 


Im 


J n ( Z ) 


K(x) 

Z 




see (2.2,25) 
see (3.2.1H) 
see (2.3.19), (2.3.22) 

amplitude of translation mode of vibration 

perturbation enthalpy 

= CTf(2x/c), airfoil thickness 
distribution 

dimensionless perturbation enthalpy 
acoustic far field enthalpy 
surface enthalpy 

/T 

unit vectors along the x, y, z axes 
Imaginary part of a complex quantity 
Bessel function 
ojc/2a, Helmholtz number 
kernel function, see (3.2.8) 
reference length 

coefficients in definition of load 
distribution, see (3.2.10) 

load coefficients for the translation 
mode, see (3.2.21) 
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L(t) 

n 

P 

P' 

P 


p ref 

p rms 

P 



Pr 


Met 

P(-1),P(1) 

M 


load coefficient for the rotation mode, 
see (3.2.22) 

surface load distribution, see (3.2.2) 
unit normal vector to C Q , see Fig. 2.1-1 
background pressure 
perturbation pressure 
unit vector, see Fig. 2.1-1 

-4 , 

dB reference pressure, 2 x 10 dynes/cnr 
root mean square pressure, see (3.2.29) 
pressure in dB, see (3-2.30 and 3*2.31) 
radiated acoustic energy, see (4.1.9) 
see (3.2.23) 

Legendre polynamial 
Prandtl number 

-1 2 

dB reference power, 10 watts 
dissipated vortical energy, see (4.1.10) 
leading and trailing edge load (experimental) 
absolute value of the quantity q 
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q* 


q + 


®mn 

R 

Re 

s 

s' 

S(ky) 

see) 

t 

T 

U n (x) 

V„ 

$ • 

S’ 


complex conjugate of the quantity q 

Fourier transform of the quantity q, 
see (4.2.1) and (4.23) 

see (3-2.23) 

perfect gas constant or Euclidean distance, see (Eq. 3-3-3) 

real part of a complex quantity 

speed ratio see (4.1.5) or arc length 

perturbation entropy 

see (3-2.5) 

see (3-2.13) 

see (3-3-10) 

time 

background temperature 

Chebyshev polynomials of the second kind, see (3-2.11) 

free stream velocity 

perturbation velocity 

dimensionless perturbation velocity 

surface work, see (4.1.8) 

h 0 /a 0 , axis of rotation 
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*o (s) 

x,y,z or Xj 

t 

a 

a o 
a, e 

r(z) 

6 

^mn 

Ae 


X 


A, y 


X 


o’ 


X 


1 


v 

P 

0 

T 

T ij 

* 


vector definition of the vibrating surface 
C Q , see Fig. 2.1-1 

Cartesian coordinates 

unit tangent vector to C Q , see Fig. 2.1-1 
ogiv/4 ( or p our i er transform variable)* 
amplitude of rotational mode of vibration 
Fourier transform variables of (x f y) 
gamma function 
os 

Kroneker delta symbol 

increment of acceleration sensitivity parameter, 

see (3.5.1) or increment of in-plane vibration amplitude, 

see (3.6.1) 

see (3.2.19) 

coefficients of viscosity 

see (4.2.4), (4.2.5) 

y/p, coefficient of kinematic viscosity 

background density 



airfoil thickness ratio, see (3.1.8) 
perturbation viscous stress, see (2.1.4) 
ir-f3, directivity angle, see Fig. 3-2-1 
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3 

Y(x) 

V 2 

Y 

J? 

(-) 

grad 

div 

curl 


frequency (rad/sec) 

dimensionless perturbation vorticity 

vorticity distribution, see (4. 3. 4) 

Laplace operator 

Cp/c y , ratio of specific heats 

see (4.4.25) and Fig. 4.4-1 for graph of 

used to denote dimensional variables in Section 3 

gradient operator 

divergence operator 

vector curl operator 
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2. ANALYSIS OF A VIBRATING SURFACE OF FINITE THICKNESS 
IN A COMPRESSIBLE VISCOUS MEDIUM AT REST 

2.1 Basic Formulation and Assumptions 

We adopt the Navier-Stokes equations as the starting point for our 
analysis. These equations are linearized around a background stationary 
fluid medium that we assume to be a perfect gas. Using primes to denote 
perturbation quantities, the linearized equations can be expressed in terms 
of enthalpy h', entropy s', and velocity vj as follows: 

J-JSi.ili.-i.fei (2.1.D 

a 2 at 8x i at W 



( 2 . 1 . 2 ) 



(2.1.3) 


where the following definitions pertain: 


av k ’ / 3v i 3 v i\ 

tJ j = X — — 6.-^ + y( — ~ * — 4- Viscous Stress Tensor 

10 ax k 1J V ax 1 ax 1 / 


(2.1.4) 



Prandtl Number (2.1.5) 


Perfect Gas Law (2.1.6) 


The perturbation pressure can be calculated in terms of h' and s’ with the 
formula 



(2.1 .7) 
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We first assume that viscous processes are more important than heat 
conduction processes for the problem under consideration. This means that 
the Prandtl number is infinite and the perturbation entropy Is zero. It is 
consistent with this assumption to neglect compressible viscous effects. 
Such terms in the momentum equation lead to the ultimate conversion of 
acoustic energy into heat many wave lengths from the vibrating surface and 
are only important in the near field for frequencies in the gega hertz 
range. With the foregoing assumptions, the perturbation equations reduce 
to the following set; 


1 3h» J 

a 2 at 


div v* 


0 


( 2 . 1 . 8 ) 


3v' 

- — + grad h' 
3t 


- v curl (curl v' ) 


(2.1.9) 


with the pressure given by 


p* « p h» 


( 2 . 1 . 10 ) 


Justification of the Model 


The foregoing equations or their inviscid counterpart with v = 0 will 
be used throughout this report except in Section ^ where we modify the 
equations to include a constant streaming velocity. We have derived these 
equations somewhat formally without reference to the specific application 
we have in mind; i.e., analysis of the Brooks experiment. The linearized 
model is justified on the grounds that all of the near field and acoustic 
pressures measured by Brooks were simple harmonic in time. Furthermore, 
the maximum velocity on the airfoil surface (see Appendik A) was in all 
cases less than 0.5 ft/sec with amplitudes less than three thousandths of 
an inch. 

The formal assumption of an infinite Prandtl number can be justified 
as follows. We use Eq. (2.1.3) to deduce the order of magnitude of entropy 
variations relative to enthalpy variations. Suppose that h* varies on a 
length scale i and the basic frequency is u. Then 


L±H 


p*Pr*u)t 


h' 
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where 


Pr 


( 2 . 1 . 11 ) 





is the Stokes number* based on the length scale of the enthalpy 
fluctuations. From the table of basic data in Appendix A, the Stokes 
number is of order 3000 based on the semi-chord. If we use the radius of 
curvature of a 12$ thick airfoil, the Stokes number is still of order 30. 
The relative contribution of the enthalpy and entropy to the pressure cah 
be estimated with (2.1.7). We have 

Jill, jh'1 - i- . 1|1 - Ih'l - (2.1.13) 

p 11 y R 1 1 oPr 


Even for a Prandtl number of unity (more typical of air), the entropy 
contribution to the pressure is of order (1/o) smaller than the enthalpy 
contribution. Thus we can safely use the isentropic model (2.1.8, 2.1.9, 
2.1.10) to analyze the length and time scales that are important in 'the 
Brooks experiment. 

2.2 The Two-Dimensional Boundary Value Problem 

Referring to Fig. 2.1-1 we consider a two-dimensional surface C 0 that 
executes rigid body simple harmonic translation and rotation about the 
origin of coordinates. The surface displacement can be expressed in the 
form 


3 - Re[(ph 0 - a 0 icxx 0 (s))e iajt ] (2.2.1) 


where p is a unit vector in the direction of the translation mode with 
amplitude h D , and a Q is the amplitude of the rotation mode. The vector 
x Q (s) defines the surface C 0 where s is the arc length along the boundary 
as indicated in Fig. 2.1-1. The unit vectors tangent and normal to C 0 are 
defined by the relations 


2 

*Aetually o is a Reynolds number based on the velocity wZ and length Z. 
However, o is the more natural viscous diffusion parameter of* the problem 
and will be referred to as the Stokes number in this report, although the 
author could not find a precedent for this terminology. 
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13 


t-*a 

ds 


(2.2.2) 

3* 

R 

cN* 

X 

X* 


(2.2.3) 

The boundary value problem for the perturbation 
given by the following formulae: 

enthalpy and 

vorticity is 

2 

V 2 h* + h’ - 0 


(2.2.4) 

a 2 



V 2 fi ' - — ft' - 0 

V 


(2.2.5) 

~~ - v - = w 2 d t (s) on 

3s 3n z 

c o 

(2.2.6) 

+ v » m 2 d n (s) on 

3n 3s n 

c o 

(2.2.7) 

where 



d n (s) = (n*p)h 0 + (t*x 0 )a 0 



d t (s) ■= (S*p)h Q - (n*x 0 )a Q 


(2.2.8) 

The perturbation vorticity Q' is related to the 
the relation 

perturbation 

velocity via 

8' - i< • curl v* 


(2.2.9) 


To solve the foregoing problem, we first introduce dimensionless 
dependent and independent variables as follows: 


(x,x 0 ,h 0 ) - | (x,x 0 ,h 0 ) 



with 

k -i£ 
2a 


h' - (22-f H « O' - « 0 

Helmholtz Number 


0 



Stokes Number Based on Semi-Chord 


oe 


iir/*» 


The dimensionless problem becomes 


with 


where 


V 2 H + k 2 H = 0 


V 2 fi - a 2 n ■= 0 


3K . 1 an _ , . 

to ^to °n <s) 


3h 1 an n . . 

3s 2 3n t 


on C D 
on C Q 


D n (s) - (n*p)h Q + (?*x Q )a 0 
D t (s) - (?*p)h Q - (n*x Q )a 0 


( 2 . 2 . 10 ) 

(2.2.M) 

( 2 . 2 . 12 ) 

(2. 2. IB) 

(2.2. U») 

(2.2.15) 

(2.2.16) 
(2.2.17) 


(2.2.18) 


Finally, we introduce a conformal map from the physical plane z = x + iy to 
an arbitrary new plane w = u + iv via the transformation 


z = F(w) 


(2.2.19) 


The problem in the w-plane can be expressed as follows 
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V 2 H + k 2 g 2 H - 0 


( 2 . 2 . 20 ) 


with 


v 2 n - o 2 g 2 n « o 


v 2 . 4 * 4 

au 2 av 2 


The boundary conditions on C^, the image of C 0 


H + = n fol „ n 

3N 2 as D N (S) ‘ g ° D n 


3H _ J_ 3Q 
as ” q 2 3N 


D t (S) - g 0 D t 


( 2 . 2 . 21 ) 

( 2 . 2 . 22 ) 

under the conformal map, are 

(2.2.23) 

( 2 . 2 . 2 * 1 ) 


The connection formulae that relate the z-plane and w-plane are 


dz 

dw 


- | F* (w) | , g Q ■= g(u 0 ,v Q ) 


(2.2.25) 


2 o " F(w o ) 


, dS 


ds 


on Cl 


(2.2.26) 


as e ° 3s 


’ 3N g ° 3n 


(2.2.27) 


2.3 The Vibrating Elliptic Airfoil 


We consider the special case of an elliptic cross section that 
vibrates along the y-axis (p ■ j) and rotates about the origin of 
coordinates (see Fig. 2.1-1). Furthermore, we construct the solution in 
the w-plane where the ellipse becomes the unit circle. The appropriate map 
is the Joukowski transformation, 


z 


F(w) 


1 +T 

_ • 

2 


w + 


1-T 


2 


1_ 

w 


(2.3.1) 
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where i is the thickness ratio. Thus 


g Q (8) “ (sin 2 0 + t 2 cos 2 6) 1/2 


( 2 . 3 . 2 ) 


= T COS 0/g- 


, n*j = sin 0/g f 


■+ ( 1-t 2 \ sin 20 ■+ -*• 

-X ° ■ j-i7~ • 


(2.3.3) 


so that 


D n - h Q sin 0 - 1 — — J a Q sin 20 


D? “ h Q t cos 0 - a Q T 


(2.3.4) 


Consider the solution of (2.2.21) in the neighborhood of the surface 
of the unit circle (C^) in the w-plane. Let 


r * 1 + R/o 


(2.3.5) 


so that 


,2 [-21 . 4- - 2 - . -i- 

3r2 0+r 3r (o+R) 2 30 2 


( 2 . 3 . 6 ) 


g 2 - g 2 (0) + 0(1 /o) 


(2.3.7) 


For o ■* ® with R = 0(1), we obtain the approximate formula 


3R^ 


( 2 . 3 . 8 ) 


with the solution 
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Thus 


a = a 0 (e)e 


-g 0 e 1 ” /l| R 


-og (r-1) 

* a n (e)e 0 


9a 

9r 


r* 1 


“go n o (e) 


and from (2.2.24) we obtain 


(2.3.9) 


( 2 . 3 . 10 ) 


«o 1_ 9H + _^T_ 

0 2 “ ag Q 90 ag Q 


With the last result, the boundary condition (2.2.23) becomes 
circle plane) 


(2.3.11) 


(in the 


Furthermore, 


we 


9H _ 1 9 / 1 9H\ n 

9r a 90 \g 0 90/ N 


use (2.3.4) to show that 


( 0 ) 




The final boundary value problem for H becomes: 

V 2 H + k 2 g 2 H * 0 



( 2 . 3 . 12 ) 


(2.3.13) 


(2.3.14) 

(2.3.15) 


We remark that the last result can be derived for an arbitrary cross 
section. However, the basic point of this section concerning the nature of 
the edge loading can be made with the special case of an eliiptic cross 
section. It is important to note, however, that the coupled problem for 
enthalpy and vorticity has been reduced to a more simple problem for the 
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enthalpy alone with a higher order boundary condition due to viscosity. 
Without the second derivative with respect to 9 in the boundary condition 
(2.3.15) the solution of (2.3.14) is in general not uni que . 

The solution of (2.3.14) can be constructed for arbitrary Helmholtz 
number k. However, we are trying to deduce the character of the edge 
loading with the viscous model. It is sufficient to illustrate the nature 
of the edge load for the compact surface for which k tends to zero. Thus, 
we consider the solution of Laplace’s equation in the circle plane with the 
generalized viscous boundary condition (2.3.15). Assume a series solution 
of the form 


H 


00 

n-1 


sin n 6 
„n 


(2.3.16) 


where the coefficients are to be determined by the surface boundary 
condition. Substitute (2.3.16) into (2.3.15), multiply by sin m 6, and 
integrate over 0 < 0 < 2ir to Obtain the following system of equations: 

N 

]C C mn A n = G m m “ 1 ’ 2 ’“* N (2.3.17) 

n=1 


The matrix coefficients and surface enthalpy for the translation and 
rotation modes are defined below. 

Translation 


C 


mn 


m 


H 0 (x) 


SB £ 

°mn 


ml 



2a' £ F 2(m-n) + F 2(m+n-1) 

+ 2a [ p2m + F 2(m-1 ) j 
A n sin (2n~1)0, cos 0 = x 


(2.3.18) 


(2.3.19) 


( 2 . 3 . 20 ) 
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Rotation 


n 


mn 


3 mn a r 2(m-n) r 2(m+n) 


( 2 , 3 * 21 ) 


(1-T 2 ) 


m 


6„i + - F, 


4 J m1 T a 2m 


( 2 . 3 . 22 ) 


N 

H°( x ) Ajj sin 2n6 » cos 6 - x 

n*1 


(2.3.23) 


The coefficients F n are defined by the following integral relation: 


?„ct) - i y 2 ” — 

" " e o 


co s n6 

(e) 


de 


f/ 

•'n 


cos n6 d0 


0 (sin 2 6 + t 2 cos 2 6) 1/2 


, n-0,2,4... (2.3.24) 


These coefficients can most easily be calculated with the following 
recursion formula 


.2 


r n+2 


2(2/k»-1)F n - (1-1/n)F n _ 2 


1 - 1/n 


(2.3.25) 


and the starting values 


F = - K(k’) 


1 


F o ' ~!t 2 e(1<,) 
° irk ,z 


( 2 . 3 . 26 ) 


where 


k* « /1-i 


(2.3.27) 
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and EOO, KCk 1 ) are the standard elliptic integrals (Ref. 7); i.e., 


E(k») 


K(k') 



/1-k' sin 0 


d9 

✓l-k^sin 2 e 


de 


- (2.3.28) 


2.4 Incompressible Viscous Load Calculations 

The matrix Eq. (2.3.17) was solved numerically for unit translational 
and rotational modes 6f vibration of an elliptic airfoil section. 
Thickness ratios, x * ,05, .1, .2 , and Stokes numbers o * 10, 100, 400 
were considered. The results for the modal surface load distributions 
(2.3.16) are presented in Figs. 2.4-1 through 2.4-4. The effect of 
thickness on the load distribution is relatively small for t < 0.2. Thus, 
we have only presented results for all three values of x for the lowest 
value of Stokes number; i.e. 0 = 10 in Figs. 2.4-1 and 2.4-2. For 0 - 10 
the out of phase (imaginary) part of the load is relatively large and there 
is a tendency toward a large edge loading. This result is in qualitative 
agreement with the flat plate viscous calculations of Ref. 5. For 0 = 100 
and 400 the out of phase load becomes much smaller, thus indicating that 
less work is done against the fluid (see Section 4.4). The more important 
point is that the edge load becomes elliptic in agreement with the result 
one would obtain by applying the Neumann edge condition to the inviscid 
calculation. We conclude that with the present method of analysis the 
combined effect of viscosity and geometric thickness does not lead to a 
flat edge load distribution of the type that Brooks measured (see Appendix 
A, Fig. A. 2). In the following section we develop the solution of the 
inviscid problem with Neumann edge condition and compare our results with 
those of the experiment. In Section 4, we reexamine the viscous problem 
for the case of zero thickness. 
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ure 2.4-lc - Translational Mode: T=0.2 , 0=10 












il Mode: r=0.1, a=100. 







Lonal Mode : T=0 . 1 , 0 = 400 







3. THREE-DIMENSIONAL INVISCID LOAD AND ACOUSTIC ANALYSIS 
OF A VIBRATING SURFACE 


3.1 Basic Formulation 


Consider the problem of a rigid three-dimensional airfoil section that 
vibrates harmonically with small amplitude in a compressible inviscid 
medium at rest. Referring to Fig. 3.1-1 » the linearized mathematical 
problem can be Stated as follows (see Eqs. (2.1.8) and (2.1.9)): 


■i^h* + div v* - 0 , h* - p'/p (3.1.D 

a 2 


imv 1 + grad h’ * 0 

where 

n|*v' - imng»3^ , |x| £ c/2 

|y| £ b/2 



and 




. * _ ? dh(x) 

± K 2 “ 

2 dx 




i<(h 0 - a Q x) + t d 0 


( 3 . 1 . 2 ) 


(3.1.3) 


(3.1 .4) 


(3.1.5) 


All perturbation quantities h* , v’ , cf| are complex. The real time 
dependent quantity, q, that corresponds to any complex amplitude can be 
calculated with the relation 

q - Re(q’e iut ) (3-1.6) 
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h 0 = Vertical Translation 
d 0 s In-plane Translation 
a 0 s Rotation 

Figure 3.1-1 Coordinate System and Geometry of a Vibrating Airfoil 


^ 'd 0 

a o 

Axis of Rotation 
at x = x o 



The three rigid body modes of vibration that we consider are denoted by h Q , 
o 0 , d Q with the following meaning: 

h Q - translation along the z-axis 

a Q - rotation about the y-axis (positive, nose up) 

d Q - translation along the x-axis 

In addition to the above equations, two auxiliary conditions are in general 
required to obtain a unique solution; i.e. , 

1) Outgoing acoustic waves for |x| “ *. 

2) Continuity of h' at the edges (Neumann condition) when needed to 
establish uniqueness. 

To proceed with the solution of the foregoing problem, we first 
introduce dimensionless variables; i.e., 


(x,y,z,h 0 ,d 0 ) - j (x,y,z,h 0 ,d 0 ) 


h(x) * c t f (2x/c) 

»• -(ffH. P -(f)? (3-1.7) 

with 

k = — Helmholtz Number 
2a 


A = b/c Aspect Ratio 

t * Max h(x)/c Thickness Ratio (3*1.8) 


Equation (3.1.2) implies that h* is the potential of imv', the complex 
acceleration, so that the velocity can be eliminated from the formulation. 
Also, we transfer the boundary condition (3.1.3) to the plane z = 0±. The 
final statement of the dimensionless problem becomes: 

V 2 H + k 2 H = 0 (3.1*9) 
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z =o± “ h ° ' “° x + d o Tf,(x) l x l * r * |y| * A 


( 3 . 1 . 10 ) 


H 



for ] x | “ ® 


(3.1.11) 


H continuous at the edges when a uniqueness condition is required 

(3-V.12) 


Finally, we observe that the problem as formulated can be split into two 
parts: 

1) Anti-symmetric Enthalpy (or Pressure) Mode 


3H 
3z i 


z=0 


h o - 


“o x 


I x 1 S 1 , | y | S A 


(3-1.13) 


wi th 


H(x,y,z) = - H(x,y ,-z) 


( 3 , 1 . 1 4 ) 


2) Symmetric Enthalpy (or Pressure) Mode 


3H 

3z 


|z=0± 


+ d 0 Tf’(x) 


with 


H(x,y,z) * H(x,y ,-z) 


(3.1.15) 


(3.1.16) 


The first problem is the result of combined translation along the z-axis 
and rotation about the y-axis while the second problem is the result of 
translation along the x-axis. We treat the two problems separately in the 
following discussion. 

3.2 Solution of the Anti-Symmetric Transverse Vibration Problem 


With the boundary condition (3.1.13). the anti- symmetric solution of 
(3*1.9) that satisfies the far field radiation condition (3.1.11) can be 
expressed as a surface distribution of dipoles oriented along the' z-axis; 
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i.e. , 


H(x) 




( 3 - 2 . 1 ) 


where 


and 


L(t) * H(x,y, 0 ) * H(x,y ,0 + ) 



( 3 - 2 . 2 ) 


(3.2.3) 


4 

The acoustic far field relative to an origin of coordinates at x 0 can be 
expressed in the following form: 


-ik( |y|+?*x 0 ) 

H f (y) - - • ( *« • y ) • S(ky) (3.2.4) 

** |y I 

where 

S(ky) - ^ d 2 f e lky *^ L(t) (3-2.5) 

y - x - x 0 , y = y/|y | (3.2.6) 


The unknown load distribution L(|) satisfies a singular integral equation 
that we derive by applying the boundary condition ( 3 . 1 . 13 ) to the dipole 
representation (3*2.1); i.e.. 


where 


with 



L(t)K(5R) 


- h Q + a 0 x 


(3-2.7) 



( 3 . 2 . 8 ) 


(x‘ 


y 2 )l /2 


(3.2.9) 
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The solution of (3.2.7) is not unique. However, we are interested in 
constructing a solution that satisfies the Neumann condition at all edges 
of the airfoil section. Thus, we assume a solution of the form 


CO 

L(x) = -— v4 2 -y 2 /r-x 2 • a Q • l^U^i (x) (3.2.10) 


where U n (x) is the Chebyshev polynomial of the second kind; i.e.. 


„ , _ . sin (n+1)e 

U_(cos e) - , x - cos 0 

n sin 6 


( 3 . 2 . 11 ) 


For convenience we assume the simple elliptic spanwise dependence of the 
load in (3.2.10). Calculation of the ohordwise pressure distribution along 
the center of the airfoil is our primary concern (see discussion in Section 
3.4) 


With the assumed form of the load distribution we now derive an 
algebraic system of equations for the unknown load coefficients L n as 
follows. Substitute (3.2.10) into (3*2.7), multiply the result by 


it 2 A 2 a. 


TfF- y 2 /l-x 2 U m _.,(x) 


( 3 . 2 . 12 ) 


and integrate over the planform to obtain the following matrix equation 
for L n . 


W 

E cl 

mn n 


g m » mas1 » 2 » • • • 


(3.2.13) 


n=1 


where 




- — ■ ‘jPdxdy /A 2 -y 2 /1-x 2 U^^x) (-h 0 + a Q x) 


x o 6 m1 " 2 6 m2 m,=1 * 2 » • • • 


(3.2.14) 
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with 


x o “ h o /a o , (defines the axis of rotation, see Fig. 3.1-1) 


(3.2.15) 


and 


6 = 0 m 4 n 

mn 


1 m = n 


Also 


■'mn 


- ^dxdy • ^A 2 -y 2 /1-x 2 U m _ 1 (x) 


(3.2.16) 


_ 2 _ 

nA 


^fd^dn /A 2 n 2 /l-£ 2 U n _.,(S) K ( x “5.y-n) 


(3.2.17) 


To evaluate C^, we introduce a double Fourier transform representation of 
the Kernel function; i.e.. 


K(x,y) 


-hf 


dadg e 


i(ax+By) 


(3.2.18) 


with 


X = (a 2 + B 2 - k 2 ) 1/2 , a 2 + B 2 2 k 2 


i(k 2 - a 2 - B 2 ) 1/2 . a 2 + B 2 S k s 


and 


j ° dadB * / da / 

* •'— ® —a 


dB 


(3.2.19) 


( 3 . 2 . 20 ) 


With the Fourier representation (3.2.18), the four-fold integral over the 
planform in (3.2.17) can be reduced to a two-fold integral over the Fourier 
space. Also, we note that C ffln is symmetric and zero unless m,n are both 
odd or both even. Thus, we can split (3*2.13) into two separate problems 
for the translation and rotation modes as follows: 
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Translation 


2N-1 

52 C mn L n h) “ x o * 6 m1 ’ m-1 ,3, , . .2N-1 (3-2.21) 

n-1 ,3 

Rotation 

2N 

5Z C mn L n a) - - i s m2 • (3.2.22) 

n-2,4 


The symmetric matrix elements can be reduced to the following set of 
integrals: 


mn 


m-n 

. m.n(-l) 2 jp^k) + lQ^Ck)] 


p mn “ I S mn cosh t) 003 * 1 t sinh 2 t dt 

J n 


/ ir/2 

(k sin t) sin t cos 2 t dt 


/' 


‘tt/2 

s mn (p) - I de 


( pA sin 8) 
pA sin 8 


2 J (p cos 8 ) J n ( p cos 8) 


p c cos 8 


(3.2.23) 


The solution of (3.2.10) for the combined translation and rotation modes is 
constructed from the solution of (3.2.21) and (3.2.22) with 




(h) 


for n odd 




(a) 


for n even 


( 3 . 2 . 2 - 4 ) 
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Both the surface enthalpy distribution and the far field can be 
calculated in terms of the load coefficients L n . The resulting expressions 
are given below: 

Surface Enthalpy 


H Q (x,y) - H(x,y ,0 + ) - - L(x,y)/2 


with 



2 „2 


/A -y 


2N 

Ljj sin n 0 
n-1 



0 ■ cos 


0 < 0 £ IT 


(3.2.25) 


( 3 . 2 . 26 ) 


Far Field Enthalpy 

Referring to Fig. 3.2-1, we evaluate the far field in the center plane 
of the airfoil (y - 0) &t a fixed distance |y| from the center of rotation 
as a function of the directivity angle ip. The result is 

-ik(|y|+x 0 cos ip) 2N 

H f ^ • tanip ^ * L n *n*i n J n (k cos ip) 

4 l y l n«1 (3.2.27) 

To conclude this section we note the following formulae that enable us 
to translate any near or far field complex enthalpy distribution into a 
magnitude and phase angle and the magnitude into an r.m.s. pressure level 
in dB: 


i* 

H - | H | e 


(3.2.28) 


where H is any complex enthalpy, | H J is the absolute value of H and 4> is 
the phase relative to the surface displacement at the leading edge. When 
we compare with data of Brooks, the phase will be shifted by ir to refer it 
to the surface acceleration at the leading edge. Finally, we note that the 
root mean square (rms) pressure is given by 


p rms 


7pk 2 |H| 


(3.2.29) 
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Let 


(3.2.30) 


P - 20 log 10 (P’ ms /P ref ) (dB) 

— li 2 

where P re f * 2 x 10 dynes/cm , the standard reference pressure for dB. 
Substitute (3*2.29) into (3*2.30) to get 


P » 20 log 10 m 200 (dB) (3*2.31) 


3*3 Solution of the Symmetric Problem (Thickness Noise) 

With the boundary condition (3*1*15), the symmetric solution of 
(3*1*9) that satisfies the far field 'hadiat ion condition (3*1*11) can be 
expressed as a surface distribution of dipoles oriented along ' the' x-axis; 
i . e . ; 


H(x) 



-ikjx-f | 
e 

m 


(3*3*1) 


where the dipole strength (determined by the boundary condition (3*1*15)) 
is proportional to the airfoil thickness distribution. Thus, the solution 
of the symmetric problem has been reduced to a double integral of known 
functions. The surface and far field enthalpy distributions are evaluated 
as follows: 

Surface Enthalpy Distribution 


It is convenient to integrate (3*3*1) by parts with respect to the 
^-coordinate, assuming that f (±1 ) - O'. The resulting expression can then 
be evaluated on the actual surface to obtain 


H (x,y) m H(x,y , xf (x) ) 


■£f r( { ) i R 


-ikR 


where 




O 2 + (y-n ) 2 + t 2 f 2 (x) * 1/2 


) 1 1 ‘ 


(3*3*2) 


(3*3*3) 
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For the NACA 0012 airfoil section used in the Brooks experiment, the 
section shape function and its derivative are given by the following 
analytic expression. 

NACA Four-Digit Series Airfoils (Bef. 8) 

f ( x ) - 1.4845 /s’ - 0.63 s - 1.758 s 2 

+ 1.4215 s 3 - 0.5075 s* 4 (3-3.4) 


f'(x) * - 0.315 - 1.758 s 

/s 


+ 2.13225 a - 1.015 s- 


where 


1 + x 


(3.3.5) 


(3-3.6) 


For the Four-Digit series airfoils the surface enthalpy distributions can 
most readily be evaluated with the following integral expression: 


H 0 f ' 

* / o 

where 

R = j"(x - cos e ) 2 + (y - n) 2 + T 2 f 2 (x) I ^ 2 (3*3.8) 


Far Field Enthalpy Distribution 


(cos 6) sin 0 de 


/ 

-A 


-ikR 


dn 


(3.3.7) 


We use the integral solution (3.3.1) directly to evaluate the acoustic 
far field relative to the same origin of coordinates that we used in the 
anti-symmetric problem (see Fig. 3.2-1). Also, we evaluate the enthalpy on 
the plane y » 0. The result is 

-ik( |y |+x Q cos ip) 

H P = - ikd rt — — • cos ip • S(k cos ip) (3-3.9) 

f 0,1 i?i 
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where 


/ i0 cos 0 

e f(cos 0) sin 0 <30 (3*3*10) 

and f(cos 0) is given by (3.3.4) for the NACA Four-Digit Series airfoil 
section. 

3*4 Comparison of Flat Plate Inviscid Theory with Experiment 

The theoretical results presented in Appendix B are calculated with 
the solution of the Anti-Symmetric Problem in Section 3*2. The surface 
pressure and phase are calculated with (3*2.25) and the far field pressure 
is calculated with (3*2.27). The results are presented in exactly the same 
format as those of Brooks (Ref. 1) with the discrete symbols denoting the 
experimental results and the solid curves denoting the calculations of 
inviscid theory. The solid curve denoted "base line (A)" in each figure is 
calculated with the vibration amplitude (oq) reported in the experiment and 
the actual geometric aspect ratio A ■ 0.666. The remaining solid curves 
differ by 2 dB increments and are convenient for reading off the 
differences between theory and experiment. 

The dashed curve in each figure is calculated for an aspect ratio A* = 
0.8 and is included for two reasons. First, we can deduce the sensitivity 
of the surface and far field pressure to changes in aspect ratio. Second, 
we suggest that the calculations based on A* can be compared directly with 
the experimental results. The reason is that we hay e used a single 

spanwise pressure mode shape of the form /l-ty/A; to simplify our 

theoretical calculations. Emphasis has been placed on the representation 
of the chordwise load distribution. In Fig. 3*4-1 we compare the 

theoretical spanwise distribution with the experimental distributions (see 
Appendix A, Fig. A. 3) • For convenience we arbitrarily set the center line 

pressure to 100 dB. It is seen that the assumed theoretical elliptic load 

based on the actual aspect ratio decays more rapidly than the experimental 
load near the spanwise edges. The rational way to determine whether the 
theory will predict the spanwise load is to include more modes in the 
assumed spanwise load distribution. This task should be carried out at 
some point. A simple and somewhat ad hoc alternative is to adjust the 
single assumed spanwise mode with the aspect ratio parameter such that the 
integrals of the experimental and theoretical mode shapes are the same. 
This leads to the value A* « 0.8 and the intermediate load curve in Fig. 
3*4-1. Also, in Fig. 3.4-2 we plot the peak surface and far field 

pressure as a function of aspect ratio for k = 1.255. The center line 

surface pressure is only increased by about 1 d"B with the increase in 

aspect ratio while the far field pressure is increased by 3 dB. It is the 
authors opinion that the calculations with the single spanwise mode based 
on A* are more representative of the solution that would be obtained with 
more assumed modes. The important point is that the modif ications only 
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affect the acoustic far field in any significant way with virtually no 
change in the centerline surface pressure. 

Several qualitative differences between theory and experiment are 
noted; 

V. Since we have built in a Neumann continuity condition, the theoretical 
edge loading is elliptic while the experimental behavior is more or 
less flat to within a few percent of the edges. 

2. The surface pressure is underpredicted across the entire chord with the 
largest discrepancy appearing in the region forward of the axis of 
rotation. 

3. With the aspect ratio A* the far field pressure amplitude is predicted 
within 1 or 2 dB. We also note that if we choose the curve of surface 
amplitude that is in best agreement with the measurements, then the 
corresponding far field amplitude is in excess of the experimental 
measurements. The inviscid calculation of Brooks also showed this 
behavior. 

*J. The phase calculations are in qualitative agreement with the 
measurements. However, the theory does not predict the far field phase 
blip near 150 deg that is present in all of the experimental data. 

A summary of the quantitative differences between inviscid theory and 
experiment is presented in Table 3» 1 • The surface pressure amplitudes 
forward of the axis of rotation (leading edge) are 5 to 7 dB smaller than 
the measurements and the levels aft of the axis of rotation (trailing edge) 
are 3 to 5 dB smaller. These dB differences correspond to amplitude ratios 
of 1.5 to 2.5 and are totally unacceptable from the standpoint of an 
aerodynamic load calculation whether or not the aspect ratio correction is 
admitted. The far field sound pressure levels based on A* are only 1 or 2 
dB smaller than the measured values, and the functional form of the 
calculated amplitude is in reasonably good agreement with the experimental 
results. With the geometric aspect ratio the far field sound pressure is 
underpredicted by 3 to 5 dB. The surface phase calculations are generally 
within 5 to 10 degrees of the measured values except for the case k » 1.255 
with amplitude 82.1 yrad where the differences on the rearward part of the 
surface are 15 to 20 deg. The far field phase differences range from 10 to 
20 degrees. 
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TABLE 3.1 SUMMARY OF DIFFERENCES BETWEEN INVISCID 
FLAT PLATE THEORY AND EXPERIMENT 


Helmholtz Number (k) 

1.255 

1.989 

ct 0 ( prad) 

82.1 

259.3 

30.8 

97.6 

Surface L.E. 

5 

7 

6 

6 

Amplitude (dB) 





T.E. 

3 

5 

4 

4 

Surface 



- 


Phase (deg) 

15-20 

5-10 

+ 10 

+ 10 

Far Field 





Amplitude (dB) a* 

0 

2 

1 

1 

Far Field 





Phase (deg) 

20 

10 

10-15 

10-15 


The main conclusion of this section is that inviscid theory based on 
the measured antisymmetric vibration mode shape cannot predict the measured 
surface pressure. The far field sound pressure is predicted with 
satisfactory accuracy when calculated with A*. In the following sections 
we investigate two plausible explanations for the large differences between 
the experimental and theoretical surface pressures. 

3.5 Effect of Acceleration on Surface Pressure Measurements 


The point of view that we adopt in this and the following section is 
that inviscid theory (with the actual aspect ratio A) is capable of 
predicting the outcome of a surface vibration experiment. We then focus on 
two possible sources of experimental error. In this section we suppose 
that the surface pressure transducers measure a component of acceleration 
in addition to pressure. Since the two signals are nearly in phase, there 
must be an enhancement of measured surface pressure level. Also, since the 
acceleration increases linearly towards the leading and trailing edges, 
there will be a tendency to "flatten" out the expected edge load 
distribution. To estimate the effect of the supposed acceleration 
contamination, we add an arbitrary component of acceleration to the surface 
pressure or we can add an arbitrary component of displacement to the 
enthalpy given by (3.2.25); i.e., 

H^ e) = H 0 + e • o D • (x-x Q ) (3.5.1) 


where e is a measure of the level of acceleration contamination. 
Alternatively, we can write the transducer signal s' as a linear 
combination of pressure and acceleration; i.e., 
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s 9 * p 9 + e s e 


ujcX? 

2/ 


/we 

P ( 


Hq + €« 


CuT 

2 


o 0 (x- x 0 ) 


CHq + e * °0 (x " x 0 )] 


(3.5.2) 


with 


e* 



e 


* 2.76 Pascals/g 


(3.5.3) 


The surface and far field calculations, based on the modified surface 
pressure (3.5.1) are compared in Appendix C. The conclusions are as 
follows. The peak surface pressure levels are clearly much closer to the 
edges (note the dashed lines), and the edge pressure distribution is much 
flatter when an acceleration component is added to the base line pressure. 
The value of e that gives the best agreement of measured and calculated 
surface pressure data is about .4 to .5 or e* * 1.08 to 1.35 Pascals/g. It 
should also be noted that the' postulated acceleration contaminant brings 
the surface signals both fore and aft of the axis of rotation into much 
better agreement with the data. The far field sound pressure levels 
corresponding to the best fit of the near field (i.e., e *= .*1 to .5) are 
greater than the measured sound by 3 to 5 dB. It is important to keep in 
mind that if the measured surface pressure 'levels were contaminated by 
acceleration then it would be necessary to remove the error before 
calculating the far field. 

In conclusion we remark that extensive calibration and acceleration 
sensitivity studies were made in preparation for the experiment (Brooks, 
private communication). Our estimated sensitivity value e* of the order of 
1 Pascal/g is considered to be much too large. 
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3.6 Effect of In-Plane Vibration on the 
Surface Pressure and Far Field 


One of the fundamental assumptions of the Brooks experiment is that 
the surface pressure amplitudes are identical on opposite sides of the 
airfoil but exactly 180 degrees out of phase. Unfortunately, the mode 
shape was measured with a single accelerometer near the leading edge that 
was oriented to measure only the transverse acceleration, the assumption 
being that the airfoil with its counterweights (see Appendix A, Fig. A. 1) 
vibrates in a pure torsional mode about the axis of the „ attached shafts. 
The maximum amplitudes of vibration are of the order of 10 - ^ inches. It 
seems possible that due to misalignment of the counterweight shafts that 
the vibration mode is actually a coupled torsion and in-plane bending 
vibration mode. This would require, of course, that the resonant frequency 
of the pure in-plane bending mode be very close to the torsional frequency. 
These speculations could be checked with a simple vibration analysis or* 
with a measurement of the in-plane acceleration. A recommendation to this 
end will be made in Section 5. Here we examine the consequences of adding 
an arbitrary component of in-plane vibration to the transverse vibration 
mode. The important point is that such a vibration mode leads to a 
component of surface pressure that is in phase on opposite sides of the 
airfoil. If an in phase component were present in the Brooks experiment, 
the corresponding surface pressure would have been treated as a lateral 
dipole because of the assumed symmetry. Thus, it is doubly important that 
both the vibration mode and the surface pressure be measured to check the 
assumed symmetry. 

To calculate the effect of in-plane vibration we combine the results 
of Sections 3.2 and 3*3. Thus, 


+ e H 0 ^ Surface (3.6.1) 

+ e H^ a) Far Field (3-6.2) 

are given by (3.2.25) and (3*3.7), respectively, and 
given by (3.2.27) and (3*3.9). The factor e in the 
is the precentage of the in-plane vibration mode. The 
results of the combined calculation are presented in Appendix D. The base 
line curve with e = 0 is indicated on each graph. The successive curves 
are presented for increments of e * 0.2. Note that a value of e - 1.0 
means that the amplitude of in-plane vibration is equal to the lateral 
displacement at the trailing edge due to the torsion mode. The 
calculations were carried out for the NACA 0012 airfoil. The pressure 
singularity at the leading edge is due to the infinite slope of the 
thickness distribution and could be removed. The significant point is that 
the surface pressure levels can be increased by in-plane vibration with the 
largest increase forward of the hinge axis. On the other hand, the far 
field amplitude is not changed appreciably except for relatively large 


H o H o 


H - 


wfcere EX §n<3 H 


l ( 4 ) or* H 

O 

are 

present calculation 




‘O 

and 


H£ a) 



values of in-plane displacement. Also, the surface phase is relatively 
insensitive to the in-plane modal amplitude but there is an effect on the 
far field phase. 

It is not intended that the calculations presented in Appendix D be 
compared directly with the experimental results to draw concrete 
conclusions. We do not know what portion of the modal amplitude, if any, 
is due to in-plane motion or what its phase may be relative to the 
torsional motion. We present these results to draw attention to the fact 
that in-plane vibration could lead to an enhancement of' the surface 
pressure with a much flatter edge loading like the one that was measured. 
A vibration analysis and further measurements, will be essential to verify 
or discount the supposition considered in this section. 
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4. ANALYSIS OF A VIBRATING SURFACE OF ZERO THICKNESS IN A 
COMPRESSIBLE VISCOUS FLUID MEDIUM IN UNIFORM MOTION 

4.1 Statement of the Problem 

In Section 2, we have shown, via a viscous analysis for finite 
thickness airfoils, that the Neumann edge condition is applicable as a 
uniqueness criterion for inviscid analysis when the Stokes number o is of 
order 100 or greater. (The Stokes number is of order 3000 in the Brooks 
experiment, see Appendix A.) In Section 3. we showed that 'the inviscid 
analysis with Neumann edge conditions cannot predict the surface pressure 
measurements of the Brooks experiment. The alternative is that either 
something has been overlooked in the viscous analysis or in the experiment. 
We further discussed two possibilities for experimental error. In the 
following, we re-examine the zero thickness problem for a compressible 
viscous medium with a low speed flow that could, for example, be maintained 
by acoustic streaming. The relationship between the surface load and shed 
vorticity is our main concern. Direct load calculations will not be made. 
In addition, we examine the energy balance in the presence of flow and 
generalize the arguments developed previously (Ref. 2) that led us to 
consider the viscous problem. 

The mathematical problem is obtained by a straightforward 
generalization of (2.1.8) and (2.1.9). We replace the local time 
derivative by the linearized convective derivative; i.e., 


at 




J- + v X 

at “ ax 


(4.1.1) 


where v w is an assumed low speed flow. With the same dimensionless 
variables used in Sections 2 and 3* we obtain the following simple harmonic 
boundary value problem: 


(2 DH 
Dx 


Div t 


(4.1.2) 


with 


+ Grad H * — -r Curl (Curl $) 
Dx 0 2 


_D_ 

Dx 



i 


(4.1.3) 


(4.1.4) 
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s 


(4.1.5) 


2v„ 


(DC 


and the boundary conditions 

£•$ * - ia Q (x-x 0 ) 

&x$ - 0 


x| S 1, |y | S A 


z - 0 ± 


(4.1.6) 


Also, we derive the following second order dimensionless energy balance by 
cross multiplying (4.1.2) and (4.1.3) with H* and (complex conjugates of 
H and V) and integrating over a large spherical volume |x| < r Q . The 
result is 


ar- 




(4.1.7) 


where 


X- 


*0' o' 


dS Re(H f V r ) 


x -r 0 »1 


^--5/ ov |S| 2 

0 • / |*l<'-o 


L 

Surface Work 

(4.1.8) 


Radiated Acoustic Energy 

(4.1.9) 


Dissipated Vortical Energy 

(4.1.10) 


All symbols have the same meaning as in Sections 2 and 3. 
4.2 Representation of the Solution in Fourier Space 


Following the analysis of Section 3» we first derive a representation 
of the solution of Eqs. (4.1.2) and (4.1.3) in Fourier space. These 
results are used in the subsequent section to derive a relation for the 
dissipation and an expression for the shed vorticity in terms of the load. 
Introduce the double Fourier transform pair; i.e., 


'-'oq * q 




oq m q(x,y) 


+ (a, 8) - j 

=1 f 


-i(ax+By) 

dxdy e q(x,y) 


i(ax+ey) + 
dadg e q (a, 8) 


(4.2.1) 

(4.2.2) 
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Then transform Eqs. (4.1.2) and (4.1.3) and seek a solution of the form 


_+ n+ ~ A I Z I 

q * q e 


(4.2.3) 


There are two modes with which we can construct a complete representation 
of the solution; i.e., 


A Q - [a 2 + g 2 + io 2 (1+as)] 1/2 Viscous 


(4.2.4) 


A, * [a 2 + 3 2 - k 2 (1+as) 2 ] 


Compressible 


(4.2.5) 


If L denotes the transform of the surface load distribution, then we have 
the following representation of the transformed solution: 

Enthalpy 


L* - X lW 

H - — — e sgn z 

2 


(4.2.6) 


Velocity 


aL + f ~V Z I - A iM) 

HuasT \ e “ e / sgnz 


(4.2.7) 


+ ( “ A 0 I Z I " A ll 2 l 


V - - K \ e 
2( 1 +os) \ 


sgn z 


(4.2.8) 


. + 2 + „2 / - A oM _A 1 1 2 1 

u + _ a p e - i o 

W 2i(1 +as) A 0 V® Al 


(4.2.9) 


Vorticity 


2 + -A rt zj 

6.V . ol 

21»o 


(4.2.10) 


ao 2 L + " A °I Z I 

• e 

2iA 0 


(4.2.11) 
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0 


(4.2. 12) 


If we invert Eq. (4.2.9) Jor the z component of velocity and evaluate the 
result on the plane z - 0 , we obtain the following: 

Integral Equation for the Transformed Load 


1 

2(2tt) 2 



i(ox+6y) 

odfi e 



«o (x “ x o ) * 

(4.2.13) 


where the boundary condition (4.1.6) has been used to evaluate the right 
hand side. 

We do not propose to solve the integral equation (4.2.13). Rather, 
our objective is to calculate the shed vorticity and dissipation in terms 
of the load distribution and to provide an energy budget for the 
experiment. 

4. 3 Distribution of Spanwise Vorticity 


In this section we derive an expression for the total spanwise 
vorticity that is shed from the trailing edge in the presence of a low 
speed flow for a load distribution of the form 


L(x,y) - A 2 - y 2 • L(x) (4.3.1) 

TT PT 

with 

L(x) - P( x ) (1-x) e (4.3.2) 


where P(x) is bounded near the trailing edge for some value of e. Thus 


P(1) 

“ Lim 
x+1 

L(x) 

for some e 

(4.3.3) 

(1-x) e 

Furthermore, we calculate 
spanwise vorticity; i.e.. 

the 

following 

weighted average of 

the total 
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Y(x) 


(4.3.4) 


where 




Y y dy 


f 

*■ 00 


Sly dZ 




O 


With (4.2.11) we get 


iao 2 L + 


(4.3.5) 


(4.3.6) 


and with (4.3.1) and (4.3.4) we obtain 


Y(x) 


3xAo 2 3 f lax 

4 3x J 


L + ( a)da 



(4.3.7) 


For large Stokes number (o » 1) the double integral can be evaluated 
asymptotically; i.e., 


Y(x) 



L*( oOdot 
a 2 +io 2 (1+as) 


(4.3.8) 


The last integral can be expressed as a convolution of the load with 
exponential functions that correspond to the roots of the denominator. The 
final result for the total vorticity downstream of the trailing edge is 


where 


Y(x) 


oA 

p.q 



-A(1-0 

L(C)e d£ 


— A( x— 1 ) 
e 


x > 1 (4.3.9) 


A = op(q-l) 


( 4 . 3 . 10 ) 
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p 


(4.3.11) 


os _ J_ Y <® 
2 2 /ojv 


q * 



(4.3.12) 


The last result has been derived under the condition that o » 1 and 
p * 0(1). Finally, we use the asymptotic form of the load distribution 
near the trailing edge (4.3.2) to evaluate (4.3-9) for large X. We get 


Y(x) = - — • P(1 ) *e 

p.q 


-X(x-l) Z* 1 

s J, 


(1-£) e e 


-A( 1-C) 




oKimin) c ~x(x-d ^ xk1 
P-q X E 


(4.3.13) 


where r(z) is the gamma function. The total vorticity downstream of the 
trailing edge is 


T 



Y(x)dx 


J_ P(1)r(H-e) 
o e P*q[p(q-1 D 1+E 


(4.3.14) 


For all edge loads that tend to zero (i.e., e > 0) the total shed vorticity 
decreases with increasing Stokes number* and is of 0(1 /o e ). For. a constant 
edge load (i.e., e * 0), the total shed vorticity is only a function of the 
speed ratio and is of 0(1). 

Typical results for the decay of wake vorticity are presented in Figs. 
4.3-1a, 1b, 1c for p =* 0,1,2. The abscissa is the product of Stokes number 
and distance downstream in semi-chords. For the Brooks experiment with 
o «= 3000, the actual downstream extent of the entire wake structure is of 
0(.002, .016, .1) semi-chords, respectively for the three values of p. If 
the assumed streaming velocity is of the order of the trailing edge 
velocity of vibration, then the experimental value of p is at most about 
0.35 based on data in Appendix A, Table A. 2. Thus, the -entire wake 
structure decays within a distance of about .01 semi-chords. 
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Figure 4,3-2 - Variation of Amplitude and Phase of Shed Vorticity; 
(Compare Figure 4. 3- lb). 



The results in Fig. 4.3-2 are presented in the form of amplitude and 
phase for p - 1 and e « 0. The amplitude decays exponentially downstream 
while the phase increases linearly. Thus the wake can be envisioned as a 
train of discrete vortices that convect at constant speed, v„, and decay 
(dissipate) with increasing distance. The magnitude of the first or 

starting vortex is proportional to o'” e » Thus, we conclude what is 
intuitively obvious that a flat edge load (e - 0), if it can exist, must 
have a vortex that is /d stronger than a trailing edge that is elliptically 
loaded (e - 1/2). 

4.4 Energy Balance 


In Ref. 2 we presented an argument for the internal consistency of the 
Brooks surface and far field data. We estimated the viscous dissipation 
and far field acoustic energy. For' the Helmholtz numbers of the Brooks 
experiment and the measured load distributions, the surface work is almost 
equally partitioned into sound and vortical energy (see Ref. 2) . The 
argument is based on two-dimensional analysis and the conclusions based on 
the energy argument led to our subsequent studies of the surface load with 
the viscous theory. As reported in Section 2, the theory has been 
unsuccessful in reproducing the Brooks edge loading and always leads to the 
same conclusion; i.e. , the Neumann condition with inviscid theory is a 
valid approach. The puzzle still remains, however, and we want to 
reexamine the energy argument in this section with the complete 
three-dimensional compressible viscous theory. In particular, we estimate 
the surface work the radiated acoustic energy and the dissipated 
vortical energy y . The appropriate dimensionless formulae are given by 
Eqs. (4.1.8), (4.1.9), and (4.1.10). The dimensional work or energy can be 
calculated with the formula 



3 2 

pa J £ * 

2 2 


non-dim 


(4.4.1) 


where,^ !110n-< ^ m is either of or ^ . Also, the power level in dB 

is given by 


(dB) 


10 log 10 


j^dim 



where J^f ef 


is the standard reference power of 


10 12 watts. 


(4.4.2) 
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Surface Work, 


form 


For the theory and experiment, we adopt a load distribution of the 


L(x,y) - A-y 2 L(x) |x| < 1 , ]y| < A (4 . i|. 3) 

ttA 2 


Substitute (4.4.3) into (4.1.8) to obtain 


asr-%/ 

• / -i 


(x-x Q ) Im L(x)dx 


(4.4.4) 


For the theoretical load distribution we have (see (3-2.6) ) 


L(x) 


“b Y, Vi 


(x) 


n=1 


so that 


f • 4 (" x o Ira L i + j Im h ) Theop y 


(4.4.5) 


(4.4.6) 


For the flat edge load of the experiment we assume 


w 

L(x) = L n P n _-,(x) 


n=1 


(4.4.7) 


where P n (x) is the Legendre polynomial. Thus, we obtain 

* 2a 2 ^-x 0 Im L-j + — Im L 2 ^ Experiment 


(4.4.8) 


where and L 2 are to be evaluated from the experimental load 

distribution; i.e., 


63 



W x - y) 


n-1 ,2 


(4.4.9) 


ImL " ' ( n ' ijf, £ 


dy 


Radiated Acoustic Energy, 

With (4.1.9) and the fact that in the far field 


V r « k H f 


we obtain 


k /,, di! r i l"fl 2 

v l x h r o 

and with (3.2.4) and (3-2.5) of Section 3 we get 


(* 


(i 


H. 


k z 

fl * 4^ * ^ 


|S| 


where 


f & lky-t 


(4 


(* 


Introduce spherical coordinates 


(x,y,z) « r(cos i}i cos <p, cos i|> sin 4>, sin i|>) ( A 


and evaluate 


TfAr J- a 


ik*n cos \|> sin <J> 

ri e 


Y'L n -f ilA-f* u n .,<0 e 

" ■'-1 


ik*£ cos cos 4> 


2ira 0 J^A) 

“I A 


(l 


.4.10) 

.4.11) 

.4.12) 

.4.13) 

.4.14) 


.4.15) 
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where 


Ev 


n*i 


■Ue) 


n n 


(4.4.16) 


n=1 


with 


A - kA cos ij> sin <j> 
0 • k cos ip cos <f> 


(4.4.17) 


The final formula for the sound power is 


/ ir ~ir/2 J (A) 2 

d<f> / cos ip sin 2 i|»di|> — |Q| 2 (4.4.18) 


Note that the sound power depends on all of the load coefficients (see 
(4.4.16)) in contrast to the surface work. 

Dissipated Vortical Energy, 

We turn now to the calculation of the energy that is dissipated in the 
form of vorticity. From (4.1.10), and the Parseval relation 


we obtain 



(4.4.19) 


& 

v 


1 

(2tt) 2 o 2 



(4.4.20) 


From (4.2.10) and (4.2.11) we have further 



oW+B 2 ) I L*| : 



e 


- 2 ReA 0 |z| 


(4.4.21) 
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so that the integration over z can be carried out in (4.4.20) to obtain 

& « — g 2 - qfcodg (a2+&2) l L+ l 2 (4.4.22) 

4(2ir; J |A 0 | 2 ReA 0 


To evaluate (4.4.22) we adopt the following generalization of the surface 
load distribution with which we can simulate the theoretical or 
experimental results: 


L(x,y) - 2 -~ v4 2 -y 2 (1+x) v (1-x) e P(x) (4.4.23) 

itA 2 


where P(x) is a polynomial and e,v are to be determined by an appropriate 
limit process at the leading and trailing edges (see Section 4.3 for 
details). The asymptotic evaluation of (4.4.22) is a tedious step that we 
outline in Appendix E. The final result is, for e * v 



4 e (r(1 +e) ) 
_2e 



where 


E (os)(jP(-1)l 2 + |P(1)| 2 ) 

(4.4.24) 




(4.4.25) 


2e 

For any 0 < e < 1/2, the dissipation is of order l/o and tends to zero 
with increasing Stokes number. We can compare 'this with the results of 
Section 4.3 where we found that the shed vorticity also has a similar 
behavior. The interesting fact is that the dissipation is practically 
independent of Stokes number for the flat edge load; i.e., 


For e * 0 


^ 2/2 
V o 

3ir*A 


( j PC-1 ) | 2 + |P(1)| 2 ) 



(4.4.26) 


where 
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du 


(4.4.27) 



The function y 0 u) is plotted in Fig. 4.4-1. We see that the dependence 
on 6 is relatively weak; 2 < 4.5 Over the full range of 6. Thus, 

the total dissipation is primarily a function of the edge loads and is 
proportional to Oq. 

In Table 4.1, we give the energy budget for the Brooks experiment. 
The total acoustic sound power is estimated with (4.4.18) with a 0 chosen to 
give the experimental sound pressure levels in the far field (see 
Figs. A. 4, A. 5, Appendix A). 


Table 4.1 Energy Budget for the Brooks Experiment 


k 

o Q (prad) 

Si (dB) 

(dB) 


Si 

A dB 

1.255 

82.1 

. 

86.9 

94.9 

91 .9 

95.5 

3 • 6 

1.255 

259.3 

98.9 

106.8 

103.9 

107.6 

3.7 

1.989 

30.8 

91.2 

94.4 

94.2 

96.1 

1.9 

1 .989 

97.6 i 

101.2 

1 

104.4 

104.2 

106.1 

1.9 


The total surface work is estimated indirectly with the inviscid acoustic 
calculation. The surface work is equal to the sound power radiated when 
the inviscid formula is used. Thus, we use (4.4.18) with a 0 chosen to give 
the far field sound pressure levels calculated by Brooks (see Figs. A. 4, 
A. 5, Appendix A). The viscous dissipation is estimated with (4.4.26) where 
the value of Jr§ corresponding to the no flow condition is used, i.e., 
^0 ■ 2.2. see Fig. 4.4-1. The difference between the surface work and the 
sum of the sound power and dissipation is 3.6 dB for k « 1.255 and 2 dB for 
k * 1.989. These results are in qualitative agreement with our initial 
estimate of the energy partitioning in Ref. 2. As a matter of interest, we 
also show in Fig. 4.4-2 the results for the surface work and sound power of 
a completely inviscid calculation over the Helmholtz number range from 
0.3 < k < 4.0 at a constant amplitude of 30.8 pradians. The energy is 
conserved (maximum error - 2 dB) and grows monotoni cally with increasing k. 
These results provide an indirect check on the spectral method of 
calculation used in Section 3. 
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Figure 4.4-2 - Surface Work and Sound Power Versus Helmholtz Number 
k for the Inviscid Calculation; a =30.8 y rad. 




In conclusion, we reiterate the fact that even though the viscous 
energy estimates can be used to argue the internal consistency of the 
surface and far field data of Brooks, all attempts to directly calculate 
the surface load have led to the conclusion that a flat edge loading is not 
possible. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


A comprehensive analytical study of the Brooks experiment with a 
vibrating airfoil has been completed. The analysis is based on linearized 
theory. The effectsof viscosity, compressibility, and finite thickness of 
the vibrating surface have been investigated. The conclusions of our study 
are summarized below. 

5. 1 Three-Dimensional Inviscid Theory with Prescribed 
Neumann Edge Conditions 


1 . With he modified aspect ratio A*, the calculated acoustic far field 
pressure amplitudes are within 1 or 2 dB of the experimental 
measurements. With the geometric aspect ratio the differences are 4to 5 
dB. 

2. The calculated phase of the far field pressure is generally within 10 
to 15 degrees of the experimental value, 

3. The calculated peak surface pressure amplitudes are under-predicted by 
3 to 7 dB, with or without the suggested aspect ratio corrections. The 
calculated pressure decays elliptically near the leading and trailing 
edges in accordance with the prescribed Neumann condition. The 
experimental pressure distributions are nearly flat near the edges. 

4. The calculated phase of the surface pressure is within 5 to 15 degrees 
of the measured value . Also, the surface pressure is nearly in phase 
with the local surface cceleration over the entire chord. 

5. Inviscid theory with prescribed Neumann edge conditions cannot predict 
the measured chordwise surface pressure distribution either in 
magnitude or in functional form. The theory can predict the form and 
phase of the far field sound pressure, but the amplitude is somewhat 
sensitive to the choice of spanwise surface pressure mode shape. With 
a suggested aspect ratio correction, the far field amplitude can be 
predicted with 1 or 2 dB. 


5.2 Viscous Theory 


1 . An energy theorem is derived that balances the work done by the 
vibrating surface against the energy converted into vorticity and 
sound. 

2, With the flat edge loading and Helmholtz numbers of the Brooks 
experiment the energy is nearly equally partitioned into sound and 
vorticity. With the elliptically loaded edge nearly all of the surface 
work is converted into sound. 
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3. A postulated small streaming velocity does not alter the last 
conclusion. However , the vorticity produced is convected downstream 
while it decays . 

1 i. All of the viscous load calculations with or without finite airfoil 
thickness hav e led to the same conclusion. For all Stokes numbers 
(o *= c/2 Aj/v j greater than approximately 100 the calculated viscous 
load distribution decays elliptically and' agrees with the inviscid 
elliptic distribution near the edges. 

5. If the viscous theory is correct, then for Stokes numbers o of order 
3000, as in the Brooks experiment, the inviscid theory with prescribed 
Neumann edge conditions is completely adequate to predict the outcome 
of a surface vibration experiment. 

Alternative * 


Either the viscous theory is inadequate to predict the measured 
surface load distribution or there is a source of error in the measured 
surface load that has not been detected. 


5.3 Possible Sources of Experimental Error 


1 . The effect of contaminating the inviscid surface pressure calculations 
with acceleration was considered. The effect of acceleration is to 
increase the magnitude of the surface pressure with very little change 
in the phase. Also, the pressure peaks calculated with inviscid theory 
are shifted closer to the edges and the edge load distribution becomes 
more flat. However, to account for the 3 to 7 dB differences between 
the measured and calculated surface pressure amplitude, an acceleration 
sensitivity of approximately 1 Pascal per g would be necessary. Based 
on the calibration studies of Brooks where the question of acceleration 
contamination was carefully considered, a sensitivity of this magnitude 
is not believed to be possible. 

2. A basic assumption in the experiment is that the pressure is 
180 degrees out of phase on opposite sides of the airfoil and that the 
vibration mode shape is a pure torsional oscillation. The effect of 
violating these conditions was investigated. The effect of adding in a 
component of in-plane vibration was calculated. The surface pressure 
amplitude is increased and the distribution becomes flatter near the 
edges. The pressure on opposite sides of the airfoil due to the 
in-plane mode of vibration is in phase. The far field sound pressure 
field is affected but to a lesser degree than the near field. To 
account for the 3 to 7 dB difference between the calculated and 
experimental surface pressures, an in phase vibration amplitude of the 
order of the trailing edge displacement due to torsional vibration 
would be necessary. 
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Recommendat ions 


1. It is recommended that the experiment be partially repeated with the 
following changes: 

a. Mount additional accelerometers on the vibrating airfoil to 
ascertain the precise vibration mode shape. In particular, 
determine whether a component of in-plane vibration is present. 

b. Measure the surface pressure on both sides of the airfoil near the 
leading and trailing edges. 

c. Reexamine the acceleration sensitivity of the surface pressure 
transducers. 

2. To support the experimental measurements, it is also recommended that a 
vibration analysis be completed to calculate the frequencies of the 
in-plane bending vibration mode and the torsion mode. If these two 
frequencies are close together, it can be expected that any 
misalignment of the torsion shafts would lead to a coupled bending 
torsion vibration mode shape even though an attempt is made to excite 
only one mode. 

3. Finally, it is recommended that the inviscid surface and far field 
calculations be repeated with more assumed spanwise pressure modes. 
The ambiguity associated with the choice of aspect ratio in Section 3. 1 * 
can then be resolved. 


73 



APPENDIX A 

BRIEF DESCRIPTION OF THE BROOKS EXPERIMENT 
AND SUMMARY OF RELEVANT DATA 

The Brooks experiment (Ref. 1) was designed to check the validity of 
the inviscid Kirchhoff integral relation 


p*(x) 



P'(?) 


3 e + 

an [Hi 


ik|x-y| 

to>p*n ( ?) - e - v;r 
|x-y| 


(A.1) 


where p'(y) and v^(y) denote pressure and velocity (complex amplitude) on a 
vibrating surface S and p'(x) is the pressure in the free field. A segment 
of a NACA 0Q12 airfoil was oscillated in an anechoic chamber at frequencies 
of 301 and 477 Hz (see Fig. A. 1). The surface pressure and velocity and 
the far field acoustic pressure were measured simultaneously. The far 
field was also calculated with (A.1) using the measured surface quantities. 
The main results of the measurement and calculation are summarized irt 
Figs. A. 2, through A. 5. The data required to perform the calculations in 
the present report are summarized below: 

Geometry (See Fig. A.1) 


AIRFOIL TYPE 

NACA 0012 

Chord (c) 

1.5 ft 

Span (b) 

1.0 ft 

Aspect Ratio (A) 

0.666 

Axis of Rotation 

7.5 in. aft of the leading edge 

Location of Transducer 
Array 

1.25 in. outboard of center span 

Radius of Microphone 
Array from Center of 
Rotation (|y|) 

7.25 ft 

Atmospheric Data 

Pressure (p) 

1.013x10 s Pa 

(1 Pa = 10 dynes/cm 2 ) 

Sound Speed (a) 

1130 ft/sec 
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Ratio of Specific Heats (Y) 1 .4 

Kinematic Viscosity (v) 1.6x10 _1 * ft 2 /sec (68°F) 

In addition to the above data, we summarize in Tables A. 1 and A. 2 the 
frequency data and some information on peak amplitudes and' trailing edge 
conditions. The data in Table A. 3 are from Ref. 1. 
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TABLE A.l Frequency Dependent Parameters 


Frequency 

Helmholtz 

Stokes 

Streaming 

Hz (Rad/sec) 

Number 

Number 

Velocity 


k= o>c/4a 

c Ju) 

9 = 2 V V 

^ (ft/sec) 

301 (1981) 

1.255 

2578 

0.55 

477 (2997) 

1.989 

3236 

0.69 


TABLE A. 2 Amplitude Data 


Frequency 

(Hz) 

Peak Angle 

of Rotation 

a (prad) 

0 

Tra 

Deflection (e) 
(in x 10^) 

iling Edge Peak Va! 
Velocity (e) 

ft/sec (e/Vcov) 

.ues 

Acceleration(e) 

e/g 

301 

82.1 

0.862 

0.136 (.247) 

7.99 


259.3 

2.72 

0.429 (.779) 

25.2 

477 

30.8 

0.323 

0.0807 (.117) 

7.51 


97.6 

1.02 

0.255 (.369) 

23.7 
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Figure A.l Experimental Model and Torsional Vibration 
System. 
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Figure A«3 Results of Spanwise Pressure Survey - Surface Pressure 
Level and Phase as a Function of Spanwise Location for 
f*447 Hz. 
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Sound Pressure Level and Phase as a Function of Directivity 
Angle g at a Measurement Radius of 2.210 m (7.25) ft) for 
Airfoil Oscillating at f^Ol Hz. 
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APPENDIX B 

COMPARISON OF FLAT PLATE INVISCID THEORY WITH 
EXPERIMENTAL RESULTS 

The theoretical calculations presented in this Appendix (Figs. B. 1 
through B.4) are based on Eq. (3*2.25) for the surface pressure and phase 
and on Eq. (3.2.27) for the far field (FF) sound pressure and phase. A 
complete discussion of these results is given in Section 3*4* The solid 
curves on each amplitude graph are calculated with the geometric aspect 
ratio (A « 0.0667). The dashed curve is calculated with the corrected 
aspect ratio (A* » 0.8) discussed in Section 3.^The experimental amplitude 
data of Ref. 1 are denoted by the solid symbols, and the phase data are 
denoted by open symbols. The Helmholtz number, k, and micro-radian 

amplitude, Oq, are noted in the title of each graph. 
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Figure B.lc Comparison of Theory and Experiment; k= 1.255, a =82.1 yrad. 
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Figure B .3b Comparison of Theory and Experiment; k® 1.989, -a =30.8 yrad. 
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Figure B. 3c Comparison of Theory and Experiment; k=l . 989, a =30.8 yrad. 
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Figure B. 4c Comparison of Theory and Experiment, k=l .989, a Q~97 . 6 yrad. 
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Figure B.4d Comparison of Theory and Experiment, k=1.989, a Q -97.6 prad. 






APPENDIX C 

EFFECT OF ACCELERATION ON THE INVISCID FLAT PLATE THEORY 


The theoretical calculations presented in this Appendix (Figs. C.1 
through C.M) are based on Eq. (3.5.1) for the surface pressure and phase 
and on Eq. (3.2.27) for the far field. A complete discussion of these 

results is given in Section 3.5. The baseline curve on each graph is 

calculated with the actual geometric aspect ratio, A * 0.666. The 

successive amplitude curves are calculated for an acceleration increment 

Ae - 0.1 (« 0.27 Pa/g). The dashed curves denote the location of the peak 
surface pressures. The experimental amplitude data of Ref. 1 are denoted 
by solid symbols, and the phase data are denoted by open Symbols. The 
Helmholtz number, k, and micro-radian amplitude, Oq, are noted in the title 
of each graph. 
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Figure C. lc Comparison of Inviscid Theory (with Acceleration Contamination) 
and Experiment; k=1.255, a =82.1 prad. 
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Figure C. Id Comparison of Inviscid Theory (with Acceleration Contamination) 
and Experiment; k® 1*255, a =82.1 yrad. 
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Figure C. 2a Comparison of Inviscid Theory (with Acceleration Contamination) 
and Experiment; k=l . 255, a =259.3 yrad* 
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Figure C. 2c Comparison of Inviscid Theory (with Acceleration Contamination) 
and Experiment; k=1.255, a =259.3 yrad. 
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Figure C*3c Comparison of Inviscid Theory (with Acceleration Contaimination) 
and Experiment; k= 1.989, ot =30.8 yrad. 






gure C. 3d Comparison of In viscid Theory (with Acceleration Contamination) 
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Transducer Location (X chord) 

Figure C. 4a Comparison of Invlscid Theory (with Acceleration Contamination) 
and Experiment; k=1.989, a =97.6 yrad. 



c 



4J 

<o 

e 

6 

cd 

4J 

c 

o 

o 

e 

o 


G *H 
O ^ 
W 0) 

*h a 
p x 
cd w 
a 

6 ^ 
O G 
O cd 




<u 

V4 

G 

DO 

*H 


114 


Transducer Location (X chord) 


Ae =0.1 



Figure C.4e Comparison of Inviscid Theory (with Acceleration Contamination) 
and Experiment; fc- 1*9 89, a =97.6 prad. 
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APPENDIX D 

EFFECT OF IN-PLANE VIBRATION ON THE INVISCID 
FLAT PLATE THEORY 

The theoretical calculation presented in this Appendix (Figs. D.1 
through D.4) are based on Eq. (3.6.1) for the surface pressure and phase 
and on Eq; (3.6.2) for the far field. A complete discussion of these 
results is given in Section 3*6. The baseline curve on each graph is 
calculated with the actual geometric aspect ratio, A * 0.666. The 

successive amplitude curves are calculated for an increment of in-plane 
vibration, Ae, with amplitude equal to 0.2 of the transverse displacement 
at the trailing edge due to torsional vibration. The dashed curve denotes 
the location of the peak surface pressure. The experimental amplitude data 
of Ref. 1 are denoted by the solid symbols, and the phase data are denoted 
by the open symbols. The Helmholtz number, k, and the micro-radian 
amplitude, Oq, are noted in the title of each graph. 
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Ang l e Be ta (deg) 

Figure D. lc Comparison of Inviscid Theory (with in- plane Vibration Com t amin a t ion) 
with Experiment ; : k=l. 255, a =82. 1 yrad. 
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Figure D.ld Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment; k=1.255, a =82.1 yrad. 
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Angle Beta (deg) 

Figure D. 2c Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment ; k=l. 255, a =259.3 jirad. 



Figure D.2d Comparison of In viscid Theory (with in-plane Vibration Contamination) 
with Experiment; k= 1.255, a =259.3 prad. 

















Figure D. 3c Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment; k= 1 . 9 89 , a =30.8 yrad. 
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Figure B. 3d Comparison of In Viscid Theory (with in-plane Vibration Contamination) 
with Experiment; k® 1.989* a =30.8 yrad. 
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Figure B. 4a Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment; k=1.989, a =97.6 yrad. 
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Transducer Location (X chord) , 

14-0CT-83 

Figure D.4b Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment; k= 1 . 9 89 , a =97.6 prad. 
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Figure B. Ac Comparison of Inviscid Theory (with in-plane Vibration Contamination) 

with Experiment; k -1 .989 , a =97.6 prad. 
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Figure D.4d Comparison of Inviscid Theory (with in-plane Vibration Contamination) 
with Experiment; k=1.989, 0 =97.6 prad. 







APPENDIX E 

THE DISSIPATION INTEGRAL 


Consider the dissipation integral (4.4.22) of the main text; i.e. , 


where 


and 
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We evaluate for o >> 1 with os - 0(1). Let 
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Then 
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First we evaluate the transformed load, 
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Now consider the integral 
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where F(v) is a bounded even function of v for 0 •*> ». Then for 0 » 1, we 
have 
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so that (E.7) can be reduced to the following single integral: 
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with 


[> + i(1+6u)j 


/2 


Now we return to (E. 11) and evaluate L + (ou). For o >> 1 and u 
the principal contribution - to the integral is from the end points. 
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For the purpose of evaluating the integral in (E.15), we have 
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We further simplify the calculation by setting v * e so that 
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Finally, we note the following results: 
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Substitute (E.15), ( E. 19), and (E.20) into (E.6) to get 
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